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ABSTRACT 



We present high quality long slit spectra along the major and minor axes out to 1.5-2 (14-22 kpc) of three bright elliptical galaxies 
(NGC1600, NGC4125, NGC7619) obtained at the Hobby-Eberly Telescope (HET). We derive stellar kinematic profiles and Lick/IDS 
indices (Hfj, Mgb, -Fejois, ^'^5270^ ^£53355^^5406). Moreover, for NGC4125 we derive gas kinematics and emission line strengths. We 
model the absorption line strengths using Simple Stellar Populations models that take into account the variation of [a/Fe] and 
derive ages, total metallicity and element abundances. Overall, we find that the three galaxies have old and [a/Fe] overabundant 
stellar populations with no significant gradients. The metallicity is supersolar at the center with a strong negative radial gradient. For 
NGC4125, several pieces of evidence point to a recent dissipational merger event. We calculate the broad band color profiles with 
the help of SSP models. All of the colors show sharp peaks at the center of the galaxies, mainly caused by the metallicity gradients, 
and agree well with the measured colors. Using the Schwarzschild's axisymmetric orbit superposition technique, we model the stellar 
kinematics to constrain the dark halos of the galaxies. We use the tight correlation between the Mgb strength and local escape velocity 
to set limits on the extent of the halos by testing different halo sizes. Logarithmic halos - cut at 60 kpc - minimize the overall scatter 
of the Mgb-Vgsc relation. Larger cutoff radii are found if the dark matter density profile is decreasing more steeply at large radii. 
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1. Introduction 

One of the most important remaining questions in contempo- 
rary astrophysics concerns the formation and evolution of stel- 
lar populations in elliptical galaxies.' One can think of two 
globally different formation scenarios: the classic monolithic- 
collapse model for the formation and evolution of early-type 
galaxies (Tinsley 1972; Larson 1975; Tantalo et al. 1996) and the 
hierarchical merging scenario (White & Rees 1978; KaufFmann 
et al. 1993). In the former, early-type galaxies build most of their 
stars during a single event in the early universe while in the lat- 
ter a substantial fraction of the stellar populations is formed from 
multiple mergers and the accretion of smaller objects over an ex- 
tended period. 

Statistically, the extremely small scatter of the observed 
color magnitude relation of elliptical galaxies or bright cluster 
members suggests that the stellar populations are surprisingly 
homogeneous (Bower et al. 1992; Kodama & Arimoto 1997). 
This uniformity of stellar populations in ellipticals is also sup- 
ported by the internal tightness of the so called Fundamental 
Plane (Dressier et al. 1987; Djorgovski & Davis 1987; Bender 
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et al. 1992; Saglia et al. 1993). These findings are consistent with 
a monolithic-collapse. Moreover, the tightness of the Mg - cr re- 
lation for massive elliptical galaxies observed in the local uni- 
verse (Bender et al. 1993; Sanchez-Blazquez et al. 2007) but 
holding up to the intermediate redshift z ~ I (Ziegler & Bender 
1997; Bender et al. 1998) requires a picture of a short, and highly 
efficient star formation process at high redshift and passive evo- 
lution since then. 



In contrast, a large number of elliptical galaxies are observed 
to have kinematically decoupled cores and were explained as 
the result of dissipational major merging events by star domi- 
nated systems (Barnes 1992; Bender et al. 1989; Mehlert et al. 
1998). Furthermore, boxy galaxies are mostly more massive, 
more radio-loud, more strong X-ray emitters, more frequently 
disturbed, and less flattened than disky early-type galaxies, sug- 
gesting that these two subclasses reflect different formation 
mechanisms and evolution history (Bender & Surma 1992; Nieto 
& Bender 1989). Recently, Kuntschner (2000); Thomas et al. 
(2005a); Collobert et al. (2006); Bemardi et al. (2006); Clemens 
et al. (2006); Rogers et al. (2010) showed that early-type galax- 
ies in low density and in high density environments might exhibit 
different formation ages, and Sanchez-Blazquez et al. (2009); 
Matkovic et al. (2009) found evidence that lower mass galaxies 
have more extended star formation histories. 
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In standard closed-box models of chemical enrichment, the 
metaUicity is a function of the yield and of how much gas 
has been locked after star formation has ceased (Tinsley 1980). 
Therefore, the metallicity strongly depends on the dynamical 
parameters - the deeper a potential is, the easier it is to re- 
tain gas in the galaxy. Then, galaxies formed in a monolithic- 
collapse should have steep radial metaUicity gradients (Larson 
1976; Thomas 1999), correlated with the local escape velocity 
- since the latter is a measure of the local depth of the poten- 
tial well. On the other hand, major mergers, which come along 
with hierarchical merging will dilute stellar population gradients 
and their connection with the escape velocity as a result of vio- 
lent relaxation (White 1980). In hierarchical merging, however, 
particularly for high-mass galaxies a significant amount of mass 
accretion can be due to minor mergers. If this deposites stellar 
populations that originally formed in small-mass halos prefer- 
entially in the outer parts of the final massive ellipticals, then 
this might enhance stellar-population gradients and one would 
again expect deviations from correlations with the local escape 
velocity, but in the opposite direction as induced by major merg- 
ers. Accordingly, the observation of stellar population gradients 
and their connection to dynamical parameters can give crucial 
insight into the formation paths of individual galaxies. 

The comparison of measured line strengths to simple stellar 
population models (Worthey 1994; Vazdekis 1999; Thomas et al. 
2003, hereafter TMB03) yields an estimate of the age, metallic- 
ity and element abundances and the galaxy formation scenar- 
ios, formation time scales and star formation histories can be 
investigated in detail from the study of the kinematic proper- 
ties and line strength features (Bender & Surma 1992; Bender 
et al. 1992; Morelh et al. 2004). Serra & Trager (2007) show 
that when two bursts are considered, the metallicity and element 
abundances derived with such a procedure are roughly luminos- 
ity weighted, while the stellar ages can be biased towards the 
ones of the younger component. 

There is a considerable number of works which focus on 
the study of kinematic profiles, line strength indices and stellar 
population parameters in early-type galaxies (Bender & Surma 
1992; Bender et al. 1994; Fisher et al. 1995a; Longhetti et al. 
1998; Tantalo et al. 1998; Trager et al. 2000b,a; Longhetti et al. 
2000; Kuntschner 2000; Terlevich & Forbes 2002; Denicolo 
et al. 2005; Li et al. 2007; Trager et al. 2008). However, the pre- 
vious measurements are concentrated within Rel2 of galaxies. 
From a comparison of the stellar parameters within RJ^i with 
these within Rel2 of some early-type galaxies sample, Davies 
et al. (1993); Trager et al. (2000b); Denicolo et al. (2005) found 
that the outer regions of the elhptical galaxies are slightly more 
metal poor than the centers, and that ages are hkely to increase 
slighdy outwards. The same trends were detected by Fisher et al. 
(1995a). But still only 1/3 of the star mass is contained within 
Rel2- So far the exact behaviour of the stellar kinematics and the 
abundance ratios at larger radii is still poorly known. 

Extended stellar kinematic profiles are also welcome to 
probe the dynamics of the outer regions of elliptical galaxies, 
where dark matter halos start to dominate the gravitational po- 
tential. Together with other kinematic tracers, like the radial ve- 
locities of planetary nebulae (Napolitano et al. 2009) and glob- 
ular clusters (Bridges et al. 2006), the density and temperature 
profiles of X-ray halos (Humphrey & Buote 2010) or gravita- 
tional lensing effects (Koopmans et al. 2009), extended stellar 
kinematics profiles have been used to measure the dark mat- 
ter halos densities and constrain their assembly epoch (Gerhard 
et al. 2001; Thomas et al. 2007b, 2009). 



A question of particular interest is the DM halo sizes. It has 
been probed statistically using gravitational lensing, and it can 
be addressed, if accurate X-ray observations sample the very 
outer regions of ellipticals. Through gravitational lensing a sta- 
tistical value for the cutoff" radius of 64^}^ kpc for the average 
cluster galaxy with a velocity dispersion of 220 km/s has been 
estimated (Halkola et al. 2007). X-ray measurements, on the 
other hand, do not seem to detect the ends of halos (Humphrey 
& Buote 2010). 

Here we attempt to improve on the leverage of stellar dy- 
namical models, that can not constrain the mass profiles ac- 
curately beyond the last kinematic data point (Thomas et al. 
2005b, 2007a), by adding the further constraint coming from the 
measured absorption line strengths. Franx & lUingworth (1990) 
found that the local color is a function of the local escape ve- 
locity. Davies et al. (1993) discovered a tight correction between 
Mg2 and the local escape velocity Vesc in 8 galaxies. More re- 
cently, Weijmans et al. (2009) also reported a tight relation be- 
tween Mg^ and in NGC3379 while Scott et al. (2009) ex- 
plored correlations between stellar population parameters and 
Vesc in the SAURON sample. 

Dynamical measurements can only constrain the density pro- 
file of the dark matter halo out to radii for which data are avail- 
able. In contrast, as the Mgb line strength is a function of the 
escape velocity at a given radius, it probes the gravitational po- 
tential directly rather than its gradient, as it is the case with dy- 
namical measurements. Therefore, since the gravitational poten- 
tial depends on the whole density distribution (in principle out 
to infinity), the Mgb- V<,sc provides a novel tool to constrain the 
density distribution of dark matter in the outer regions of galax- 
ies. 

For this work we obtained deep long slit spectra of three 
giant local galaxies, NGC1600, NGC4125 and NGC7619 (see 
Table 1). They were chosen as representative of the class mas- 
sive ellipticals, that could be easily fit into the service mode 
schedule (Shetrone et al. 2007) of the Hobby-Eberly Telescope 
(HET) (Barnes & Ramsey 1998) at McDonald Observatory. 
These galaxies have been well studied in several works (Bender 
et al. 1994; Fisher et al. 1995a). The galaxy NGC1600, has boxy 
isophotes (Rembold et al. 2002) and does not show rotation both 
along major and minor axes. The galaxy NGC4125, has dust 
and ionized gas aligned along the stellar major axis (Kim 1989); 
Goudfrooij et al. (1994) also found that the dust lane is inclined 
by 10°with respect to the apparent major axis; in addition, this 
galaxy shows strong soft (0.1-2.4 keV) X-ray and far-infrared 
{A = 60 iim) emission (Condon et al. 1998). For NGC7619 
Rembold et al. (2002) found evidence of tidal interactions with 
the galaxy NGC7626. Recently, Fukazawa et al. (2006) analyzed 
Chandra data and obtained that NGC1600 and NGC7619 show 
a positive temperature gradient towards the outer parts. On the 
other hand NGC4125 has a dec fining temperature profile. Our 
spectra extend out to 1 .5 Rg for NGC1600 and more than IRg for 
NGC4125 and NGC7619. This dataset allows us to constrain the 
general characteristics of the stellar kinematics and populations 
extending to large radii. By combining stellar dynamical model- 
ing with the measured fine strength indices we derive constraints 
on the extent of the dark matter halos in these galaxies. 

The paper is organized as follows. In Sect. 2 we describe 
the observations (Sect. 2.1) and the data reduction (Sect. 2.2). In 
Sect. 3 we present the kinematics (Sect. 3.1) and the line strength 
measurements (Sect. 3.2). Conraients on individual galaxies are 
given in Sect. 3.3. We analyze the Lick indices and derive ages, 
metallicities, and a/Fe ratios in Sect. 4. Section 4.1 describes 
the models and the method used; Sect. 4.2 presents the related 
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results; Sect. 4.3 derives colors and mass-to-light ratios. In Sect. 
5 we produce dynamical models of the galaxies. We discuss the 
adopted method in Sect. 5.1, the resulting models (Sect. 5.2), 
the constraints coming from the correlation between Mgb and 
local escape velocity (Sect. 5.3) and the resulting degeneracy 
between dark halo cut-off radii and outer density slope (Sect. 
5.4). A discussion of the results and a summary of this work are 
presented in Sect. 6 

2. Observations and data reduction 

2.1. Observations 

Long-slit spectra of NGC1600, NGC4125 and NGC7619 were 
collected along the major and minor axis in the period 2006- 
2008 using the HET in service mode and the Low-Resolution 
Spectrograph (LRS) with the E2 grism (Hill et al. 1998). For 
NGC1600 and NGC7619 the slit was centered on the center of 
the galaxies. For NGC4125, the center of the galaxy was moved 
towards one end of the slit, to be able to probe its outer region de- 
spite of its larger half-luminosity radius. The wavelength range 
covered was 4790 to 5850 A with a pixel size of 0.73A. The 
slit width was 3 arcsec, giving an instrumental broadening of 
cr,„j( = 120 km/s. Several 900 sec exposures were taken at each 
galaxy position. Moreover, 900 sec exposures of blank sky re- 
gions were taken at regular intervals. The seeing ranged from 1 .2 
to 2.6 arcsec. The resulting summed spectra probe regions out to 
nearly 1 .5 for NGC 1 600 and more than 2Re for NGC4 1 25 and 
NGC7619. The Ust of the spectroscopic observations is given in 
Table 2. 

In addition, calibration frames (biases, dome flats and the 
Ne and Cd lamps) were taken. Finally, a number of standard 
template stars (K- and G-giants) from the list of Worthey et al. 
(1994) were observed to calibrate our absorption line strengths 
against the Lick system. 



Table 2. Log of spectroscopic observation, MJ = major axis; MN = 
minor axis; SKY= sky spectrum. 



Date 


Objects 


Position 


Seeing(F^ 


2006 Sep 19 


NGC 7619 


MN1,2,SKY1 


1.26 


2006 Sep 20 




MJl 


1.97 


2006 Sep 27 




MN3,4,SKY2 


1.66 






MJ2,3 


1.34 


2006 Sep 28 




MJ4,5 


1.18 






MN4,5 


1.42 


2006 Oct 29 




MN6,7 


1.99 


2006 Nov 17 




MJ6,7 


2.17 


2006 Nov 22 




MN8,9 


1.47 


2007 Jan 28 


NGC4125 


MJ1,2,SKY1 


1.64 






MN1,2 


1.36 


2007 Feb 20 




MJ3,4,SKY2 


2.31 


2007 Feb 21 




MN3,4,5,6 SKY3 


2.54 


2007 Feb 22 




MJ5,6 


1.62 


2007 Mar 14 




MJ7,8 


2.26 


2007 May 05 




MJ9,10,SKY4 


2.21 


2007 May 10 




MN7,8,SKY5 


1.63 


2007 May 11 




MJ11,12,SKY6 


1.70 


2007 Oct 09 


NGC1600 


MJ1,2,SKY1 


1.76 


2007 Oct 16 




MN1,2,SKY2 


1.74 


2007 Oct 20 




MJ3,4,SKY3 


1.42 


2007 Oct 21 




MN3,4,SKY4 


2.15 


2007 Nov 07 




MJ5,6,SKY5 


1.66 


2007 Nov 19 




MN5,6,SKY6 


1.61 


2007 Dec 15 




MN7,8,SKY7 


1.92 


2007 Dec 16 




MN9,10,SKY8 


2.20 


2008 Jan 10 




MJ7,8,SKY9 


2.29 


2008 Feb 04 




MN11,12,SKY10 


2.57 


2008 Feb 09 




MJ9,10,SKY11 


1.49 



" The exposure time for each slit is 900 seconds. 
Position angles of our galaxies. NGC1600: [MJ, MN] are [9, 99]; 
NGC4125: [MJ, MN] are [82, 172]; NGC7619: [MJ, MN] are [36, 
126]. 



2.2. Data reduction 

The data reduction used the IVIIDAS package provided by ESO. 
The pre-process of data reduction was done following Bender 
et al. (1994). The raw spectra were bias subtracted, and divided 
by the flatfields. The cosmic rays were removed with a k - cr 
clipping procedure. The wavelength cahbration was performed 
using 9 to 1 1 strong Ne and Cd emission lines and a third order 
polynomial. The achieved accuracy of the wavelength calibra- 
tion is always better than 0.6 A (rms). The science spectra were 
rebinned to a logarithmic wavelength scale. 

As discussed in Saglia et al. (2010, Fig. 2), we also need to 
correct the anamorphic distortion of the LRS. Using stars that 
were trailed several times across the slit, we mapped their traces 
orthogonal to the wavelength direction and corrected all science 
frames by means of sub-pixel shifting using the mapped distor- 
tion. We estimate that after the correction the maximal residual 
distortion is of the order of 1.2 arcsec. 

The step of sky subtraction required particular care to min- 
imize systematic effects on the measured kinematics and Une 
strengths in the outer regions of our galaxies. We have to face the 
following problems: The HET queue scheduling operations and 
observing conditions did not always allow to acquire blank sky 
spectra after the science observations (see Table 2). Moreover 
variability in both the atmospheric absorption (due to difl'er- 
ent air masses and/or non-photometric conditions) and sky level 
were present. Finally, not always a perfectly flat slit illumina- 
tion was achieved. In order to account for these problems, we 



selected the spectra for each galaxy where a sky spectrum with 
a uniform slit illumination was available, and almost photomet- 
ric conditions were achieved, yielding the largest galaxy counts 
per pixel. To correct for the inhomogeneous slit illumination, we 
produced a 4th to 6th order polynomial model of the sky spectra 
for each column in spatial direction and subtracted it from the 
selected galaxy frames, obtaining three reference science frames 
Gr - one per galaxy. In the next step, we calibrated the atmo- 
spheric absorption and sky level of the remaining galaxy frames 
G, using the following procedure. We computed the fractional 
residuals between the scaled and the reference sUt profiles 



f? * (Gr) 



(1) 



and minimized it (see below). Here is the scaling factor of the 
noise-free (i.e. the polynomial model) sky frame taken after the 
galaxy frame G„ when available, or the average of the most uni- 
form sky frames when not. The symbol () indicates the average 
in the wavelength direction and R(r) is a function of the position 
r along the slit. Moreover, is a scaling factor that takes into 
account the different atmospheric transmissions. We determined 
fP and ff iteratively such to minimize R, which in an ideal sit- 
uation should be zero at every radii. FinaUy, we computed the 
resulting total galaxy frame G,oi as: 



G,ot = G,+UGi-f,'*SKYi)/f[^. 



(2) 



4 



S. B. Pu et al.: The outer halos of NGC1600, NGC4125 and NGC7619 



Table 1. List of the morphological classifications, B band luminosity, luminosity distance, receding velocity and effective radius of NGC1600, 
NGC4125 and NGC7619. They are taken from Faber et al. (1989); O'SuUivan et al. (2001) and NED(NASA/1PAC Extragalactic Database). 



Name 


Ra 


Dec 


Type 




Dis 


Vel(km/s) 


Rei arcsec) 


NGC1600 


04h31m39s 


-OSdllmlOs 


E3 


11.03 


64Mpc 


4681 


45 


NGC4125 


12h08m06s 


65dl0m26s 


E6 


10.80 


20Mpc 


1356 


58 


NGC7619 


23h20ml4s 


8dl2m22s 


E2 


10.58 


46Mpc 


3760 


32 



In practice, due to the non-uniformity of the sUt illumination 
function the function R(r) is not always zero, but through the 
summing in Eq. 2 the differences should average out. We can test 
the quality of the calibration by comparing the profile {Gtot) with 
available broad band photometry. Figure 1 shows some examples 
of the caUbration procedure for the spectra taken along the major 
axis of NGC7619. The different type lines indicate different slit 
spectra. The notes "NO" and "YES" indicate spectra before and 
after calibration. The slit marked with "Ref is reference spec- 
trum. The soUd and open dots show the I broad band photometry 
profile (see Sect. 5.2) along major and minor axis respectively. 
The top panel shows the ratios R of MJ3 to MJl (solid line) and 
MJ7 (dotted line). As it can be seen from the figure, the MJl and 
MJ3 matched well with the photometry profiles from the center 
to the outer parts. In contrast, the MJ7 profile has larger devia- 
tions from the photometry profile in the outer regions since the 
sky was over-subtracted on the left side and underestimated on 
the right side. The summed spectrum agrees well with the broad 
band photometry out to r » 70 arcsec. The kinematics and fine 
strength profiles derived in Sect. 3 extend out to radii where the 
agreement with broad band photometry is always good. Figure 2 
shows the comparison between the summed shifted slit profiles 
and the broad-band photometry for the three galaxies. 

The procedure described above removes the diflferences in 
the sky continuum, but does not necessarily calibrate correctly 
the sky emission lines. Therefore, following Saglia et al. (2010) 
we constructed a separate zero continuum sky emission line 
spectrum that was scaled and subtracted from the summed spec- 
tra. As a final step, the continuum was removed by fitting and 
dividing by a fourth order polynomial along the wavelength di- 
rection. 

Finally, in order to gauge the residual systematic errors, we 
repeated the extraction of the kinematics increasing or decreas- 
ing the scale factor for the sky level by the value of the rms 
difference in the normalized slit illumination function (see up- 
per plot of Fig. 1). The effect on the extracted kinematics is most 
significant in the case of NGC 1 600 and is indicated by the dotted 
errors in Fig. 3. 

3. Kinematics and line indices 

3.1. Kinematic profiles 

The line-of-sight velocity distributions (LOSVDs) and kine- 
matic parameters were determined from the continuum-removed 
spectra, rebinned radially to obtain almost constant signal- 
to-noise ratio, using the Fourier Correlation Quotient (FCQ) 
method (Bender 1990) with the implementation described in 
Saglia et al. (201 0) that allows for the presence of emission lines. 
In order to minimize the mismatching, we use the templates 
from a subset of the stellar spectra library of Vazdekis (1999). 
This hbrary contains about a thousand synthetic single- stellar- 
population spectra covering the wavelength range from 4800 to 
5470 A with the resolution of 1.8 A. We used the library with 
the age of 1.00 to 17.78 Gyr and the metallicities from -1.68 to 




-100 



-50 



50 



R [arcsec] 



Fig. 1. Examples of the sky calibration for the major axis spectra of 
NGC7619. In the bottom of the figure, the X coordinate is the distance 
from the galaxy center in arcsec and the Y coordinate is -2.5 * lg{CN) + 
25.7 (where CN is counts number). The different line types indicate dif- 
ferent spectra. The notes "NO" and "YES" indicate the spectra before 
and after the calibration of the sky respectively. The "Ref means refer- 
ence galaxy, MJ SUM is the final spectrum frame which is the sum of 
the single frames. The solid and open dots stand for the 1 band photom- 
etry profile (see Sect. 5.2) along the major and the minor axis, respec- 
tively. In the top panel of the figure, the ratios R of MJ3 to MJl (solid 
line) and MJ7 (dotted line) are presented. 



0.2. We first set all of the library spectra to the resolution of our 
galaxy spectra and find the best fitting template for each radial 
bin according to the lowest RMS value of the residual (reaching 
typically 1% of the initial flux). If emission lines are detected, 
Gaussians are fitted to the residuals above the best-fit template 
and subtracted from the galaxy spectrum to derive cleaned spec- 
tra. The kinematic fit is then redone using these cleaned spectra. 
We detect emission only in the spectra of NGC4125. 

The stellar kinematic profiles of both the major and the mi- 
nor axes of NGC 1600, NGC4125 and NGC7619 are shown in 
Fig. 3, where we show the rotational velocity, the velocity disper- 
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Table 5. The linear transformation between indices measured by us and 
in LICK systems. The second column displays the best linear relations, 
and the last column show the ratio of LICK systems to our measure- 
ments. 



Index 



Linear relation 

(Lick=axHET+b) 



RMS Lick/HET 



HP 


1.091(± 


0.096)-0.327(± 


0.134) 


0.085 


0.780 ± 


0.117 


Mgh 


1.137(± 


0.025)-0.549(± 


0.129) 


0.073 


0.967 ± 


0.0245 


Fe5015 


0.892(± 


0.105)+0.980(± 


0.663) 


0.147 


1.052 ± 


0.0276 


FeSnO 


0.879(± 


0.135)+0.499(± 


0.471) 


0.143 


1.037 ± 


0.0291 


Fe5335 


0.903(± 0.085)+0.255(± 0.321) 


0.117 


0.990 ± 0.0249 


Fe54()6 


0.974(± 


().127)+().17U± 


0.265) 


0.142 


1.066 ± 


0.0373 



R [arcsec] 
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Fig. 2. The comparison between the broad band surface brightness pro- 
files (dots) and the ones derived from our summed spectra (lines) for 
NGC1600 and NGC4125. The full line and the filled dots are profiles 
along the major axes. The dashed line and the open dots show the pro- 
files along the minor axes. 



sion and the Gauss-Hermite parameters and H/\. The curves 
are folded with respect to the center of the galaxies, filled and 
open symbols stand for difl'erent sides^ of the galaxies. The 
fl(, - Re ■ e"'^^ and - Re ■ f'^^ are labeled, where the R,, 
is the effective radius and the e is the apparent axial ratio. The 
long dashed lines show the instrument dispersion. The solid lines 
present the kinematics from the dynamic models (see Sect. 5.2). 
Table 3 gives format examples of the measured stellar kinematics 
with statistical and systematic errors. The fuU listing is available 
electronically. 

Also the previous studies of kinematics are shown in the 
plot. The circles present the data of this work, the squares show 
the previous work of Bender et al. (1994) for NGC1600 and 
NGC4125, Fisher et al. (1995b) for NGC7619. As it can be seen 
from the figures, the agreement is generally good except for the 
major axis data of NGC1600. In all cases our profiles extend 
to larger radii and show less scatter than those of previous in- 
vestigation. NGC4125 has dust and ionized gas aligned along 
the major axis (Bertola et al. 1984; Kim 1989; Wiklind et al. 
1995). We also measured the kinematics of ionized gas by fit- 

' For NGC4125 and 7619: The filled symbols show the kinematic 
profiles on south-west (SW) (receding) side along the major axis 
and north-west (NW) (approaching) side along the minor axis. For 
NGC1600: The open symbols show the kinematic profiles on the north- 
east (NE) approaching side along the major axis and south-east (SE) 
side along the minor axis. 



ting Gaussians to the emission lines of Up, [0\\\]AA 4958,5007 
and [NIJ/l/l 5197,5200. The kinematics of the gas are shown in 
Fig. 6. Table 4 gives format examples of the measured gas kine- 
matics. The full fisting is available electronically. 

3.2. Line strength indices 

Line-strength indices were defined by Burstein et al. (1984); 
Worthey et al. (1994) and redefined by Trager et al. (1998). In 
this work, the index windows follow the definition of Trager 
et al. (1998). For all of our galaxies, six line strength indices 
have been measured from the cleaned spectra along the major 
and minor axis. Before measuring the indices, our spectra were 
degraded to the resolution of Lick/IDS systems. We then cor- 
rected the indices for the velocity dispersion using template stars 
and the value for cr derived in the previous section. Finally, the 
observational data need to be corrected to the Lick/IDS system. 
To do this we observed 5 stars from the Lick/IDS Ubrary using 
the same instrumental configuration used for the science objects 
and derived the offsets between our data and the Lick/IDS sys- 
tem. The comparison of our data with the Lick standard systems 
can be found in Fig. 4. The dashed fine presents the 1:1 rela- 
tion, the solid lines indicate the Unear fit the data. The data are 
in good agreement with the Lick/IDS systems, as found in Saglia 
et al. (2010) using a larger set of Lick standards observed with 
LRS and HET at a better resolution. Only for the Hfi there is 
a small, marginally significant ofl'set, possibly related to a de- 
fect column in the CCD (see discussion in Saglia et al. 2010). 
The linear relation between indices in LICK/IDS systems and 
our measurements are listed in Table 5. We conclude that the 
deviation between our measurements and the Lick system can 
be ignored, but we take into account the RMS of the calibration 
lines into the final error budget, by adding it in quadrature to the 
statistical error of each index. 

As in Saglia et al. (2010), we do not consider the molec- 
ular indices Mg\ and Mg2 that are affect by inaccurate spec- 
tral flatfielding towards the end of the slit, where vignetting be- 
comes important. Figure 5 shows the lines strength indices pro- 
files along the major (left) and minor (right) axis of our three 
galaxies. The name and position (major axis and minor axis) 
are labeled in the figure, the asterisks show the measured line 
strengths as the function of the radii in arcseconds along the ma- 
jor and minor axes, the solid lines show the model predictions 
(see Sect. 4). Finally, as for the kinematics, we indicate the sys- 
tematic variations due to sky subtraction with dotted errors. The 
blue symbols show values from the literature (see below). The 
agree well with our measurements. Table 6 gives format exam- 
ples of the measured Lick indices with statistical and systematic 
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Fig. 3. Kinematic parameters along the major (left) and the minor (right) axis of our galaxies, the corresponding names are labeled. From bottom 
to top in each panel we show: (1) Rotation velocity, (2) Velocity dispersion, (3) and (4) Gauss-Hermite parameters and H^. The solid lines 
present the kinematic moments from dynamic models (see Sect. 5.2). The doted Unes show the systematic uncertainty due to sky subtraction. The 
curves are folded with respect to the nucleus of the galaxies, filled and open symbols stand for different sides of the galaxies. The squares show 
the data published in Bender et al. (1994). 



Table 3. Format example of the measured stellar kinematics as a function of distance from the center (positive: east, negative: west) for the difl^erent 
position angles. The full table is available electronically. 



Galaxy 


R 


PA 


Vstars =t dVstat =t 




^ Stars 


H3 






(") 


(deg) 


(km/s) 




(km/s) 






NGC1600 


-69.67 


9 


-3.8 ± 12.6 ± 


24.9 


324.4 ± 11.3 ± 17.9 


0.09 ±0.04 ±0.11 


-0.10 ±0.04 ±0.15 



Table 4. Format example of the measured gas kinematics for NGC4125 as a function of distance from the center (positive: east, negative: west) 
for the different position angles. The full table is available electronically. 



Galaxy 


R 


PA 






H/J 


fOIlIl/H/3 






(") 


(tlcg) 


(km/s) 


(km/s) 


A 






NGC4125 


-81.3 


82 


-3.4 ± 25.5 


180.9 ± 25.6 


0.06 ± 0.02 


6.0 ± 1.6 


1.7 ±0.03 



errors as a function of distance and position angle, respectively. 
The fuU listing is available electronically. 

3.3. Comments on individual galaxies 

NGC1600: This is the most distant galaxy of the sample and it is 
classified as E3. Our kinematic data extend to 70 arcsec from the 
center, almost about four times the radius of Bender et al. ( 1 994). 
The kinematic profiles show no rotation along both axes and 
high velocity dispersion. The velocity dispersion shows a steep 



gradient inside 15 arcsec, becomes flat along the major axis and 
even increases along the minor axis outwards. The galaxy shows 
a weakly asymmetric profile of Mgb, Fe5015 and Fe5335 in the 
outer part. Moreover, it has low Hj8 line strength, particularly in 
the outer parts. Because accurate measurements of HyS always 
suflfer from the contamination of emission lines or sky lines, we 
carefully checked the specfra, and no emission lines or sky lines 
were detected. Our instrumental resolution could be too low to 
detect the weak emission lines embed in the bandpass or some 
weak absorption features of other elements, such as CrI /14885 
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Fig. 3. Continued. The stellar kinematics of NGC4125. The squares show the data published in Bender et al. (1994). The solid lines present the 
kinematics from dynamic models (see Sect. 5.2). 



Table 6. Format example of the measured Lick indices as a function of distance from the center (positive: east, negative: west) for the different 

position angles. The full table is available electronically. 



Galaxy 


R 


PA 




Mgb 


Fe5015 


Fe5270 


Fe5335 


Fe5406 




(") 


(deg) 


(A) 


(A) 


(A) 


(A) 


(A) 


(A) 


NGC1600 


-69.9 


9 


0.5 ±0.1 ±0.3 


2.1 ±0.1 ±2 


4.9 ± 0.3 ± 0.3 


2.0 ± 0.2 ± 1.2 


2.7 ± 0.3 ± 1.5 


-0.4 ±0.2 ±3.1 



A, Fel /14891 A, set in the pseudo-continuum of ^8 definition 
(Kom et al. 2005; Puzia et al. 2005). 

NGC4125: This galaxy shows strong ionized gas emission 
along the major axis. The stellar kinematics were measured out 
to 1.7 ae on SW side along the major axis and 3.5 bg on NW side 
along the minor axis. As Fig. 3 shows, the kinematic profiles 
are fairly symmetric, and show a strong velocity gradient along 
the major axis; inside 50 arcsec the rotational velocity increases 
rapidly from to 150 km/s, remains remains nearly flat out to 
100 arcsec, and dechnes at larger radii. 

Furthermore, it is clear from Fig. 3 that the NGC4125 shows 
rotation along the minor axis which could indicate triaxiahty 
(Bertola et al. 1984), or be the sign of unsettled material left 
over from a past merger and visible on images of the galaxy. 
The velocity dispersion profile of the stars shows a depression in 
the irmer 8 arcsec, on both the major and the minor axes. This 
depression suggests that there is a colder substructure in the in- 
nermost center. The galaxy presents positive gradients of H/J in- 
dex both along major axis and minor axis in the center, while it 



remains flat towards outer part. The HyS profile along major axis 
exhibits weak asymmetry. 

Compared to the kinematic profiles of the stars, the veloc- 
ity curves of the gas along the axes shows stronger gradients, 
see left panel of Fig. 6, and presents a maximum of about 
240km/s at 10 arcsec. Warm gas has been found already by 
Bertola et al. (1983). Our folded gas velocity curve is not per- 
fectly axisymmetric, but less distorted than what reported there. 
Finally, the rotation curve implied by the steUar dynamical mod- 
eling is much higher (*; 300 km/s) than the measured velocities. 
Therefore the gas is not following simple circular motions and/or 
is not settled in a regular disk. The gas velocity dispersion is high 
in the inner 20 arcsec (« 160 km/s after correction for instru- 
mental broadening) and unconstrained at larger radii, where it 
is less or comparable to our spectral resolution. Figure 6 shows 
also the equivalent widths (in A) of the emission line Hj6 and 
file ratios [NI]^^5 197,5200/^8 and [OIII]A5 007/1^ in logarith- 
mic units. The data along the major and minor axes are shown 
separately. The [Nl]ii5197,5200/HyS versus [Olll].15007/H/3 di- 
agnostic diagram are shown at the bottom of right panel. As it 
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Fig. 3. Continued. The stellar kinematics along the major and the minor axes of NGC7619 is shown; the squares present the measurement of 
Fisher et al. (1995a). The solid lines present the kinematics of the dynamic models (see details in Sect. 5.2). The long dashed Unes indicate the 
instrumental dispersion. 



can be seen from the plot, the ratios [Olll]/I5007/Hy6 span the 
range from to 0.5 and the ratios [NIJ/l/15 197,5200/1^ go from 
-0.5 to 0.0. According to Sarzi et al. (2010), this is the region 
driven by LINER-like emission. 

NGC7619: Our kinematic profiles for NGC7619 extend to 
almost 1.6 a<, along the major axis and 2 b<, along the minor 
axis, further out than previous works, for example Fisher et al. 
(1995a); Longhetti et al. (1998) extended only up to 2/3 R^. Our 
data (circles in Fig. 2) agree, within the error, with the measure- 
ment of Fisher et al. (1995a) (squares). From the kinematic pro- 
files, we see that the velocity dispersion of the galaxy has a shal- 
low gradient inside 10 arcsec and remains nearly flat towards 
large radii. The rotational velocity also remains almost constant 
from 8 arcsec out to 60 arcsec. Furthermore, slow rotation along 
the minor axis reveals that the galaxy is possibly triaxial (Bertola 
et al. 1984). Concerning the lines strength profiles, the galaxy 
shows constant Hy6 index from the center to the outer regions 
and it has a steep Mg^ gradient. 

4. Stellar populations 

4.1. Model and method 

In this section, we use the stellar population models of TMB03 to 
derive the age, total metalhcity and element abundance gradients 



along the major and minor axes from the measured lines indices 
of our galaxies. The TMB03 models cover ages between 1 and 
15 Gyr, metallicities between 1/200 and 3.5 solar. Furthermore, 
the models take into account the effects on the Lick indices by 
the variation of element abundance, hence, give Lick indices of 
simple stellar populations not only as the function of age and 
metallicity, but also as the function of the ajFe ratio. The age, 
total metalhcity and element abundance can be derived from a 
comparison of selected line strength indices with SSP models 
TMB03. 

The traditional and effective method of studying stellar pop- 
ulation properties uses diagrams of different pairs of Lick indices 
(Thomas et al. 2005a). The Hfi versus [MgFeJ' pair diagram is 
selected as the best age indicator because Hfi is sensitive to warm 
turnoff stars and [MgFe] ' index is considered as the best detector 
of metallicity since it does not depend on abundance ratio vari- 
ations. Using Hp versus [MgFe]' can break the age-metallicity 
degeneracy. The < Fe >- {FeSllQ + FeSJiJii)!! versus Mgi> 
pair usually is considered as the best indicator of the abundance 
of populations. However, in this study, following Saglia et al. 
(2010) we use a simple and perhaps more accurate method- 
minimization: 

^ _ ^ {indexiiob] - indeXi[mod\)^ 

index[i\ 
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Fig. 5. The lines-strength indices along the major (left) and minor axis (right). The name of the indices are labeled and the galaxies' names are also 
noted. The dotted lines show the size of the systematic errors due to sky subtraction. The solid lines show the models (TMB03) predicted lines 
strength profiles along the axes. The blue open triangles are measurement of Sanchez-Blazquez et al. (2007). 
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Fig. 5. Continued. 
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Fig. 5. Continued. The blue open triangles showing the data taken from Fisher et al. (1995a). 
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Fig. 6. The gas kinematics both along major and minor axes of NGC4125 are shown on the left; open and filled symbols show different sides of 
the galaxy. The emission lines (in logarithmic units) along the major and minor axes of NGC4125 are show on the right. From the first row to the 
third row, we show the H/3 emission EW, the [OIII]/i5007/Hj8 and the [NI]/l/15197,5200/H/3 ratios (in logarithmic units) as a function of distance 
from the center. In the bottom panels the diagnostic plots [NI]/l/l5197,5200/H/3 against [OIII]/i5007/Hj8 in logarithmic units are shown. 
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Fig. 4. The comparisons between the values for the lines strength in- 
dices derived in this work and the standard Lick/IDS indices. The 
dashed line presents the 1:1 relation, the solid lines indicate the linear 
relation fit to the data. 



where indexiyoh] and indexi[,„od] represent the ; observational in- 
dices and model indices respectively, o-i[ob\ is the observational 
uncertainty of indices. The best fitting age, metaUicity and 
alFe can be derived by finding the minimum of all selected 
lines indices to the SSP models. We chose HjS, Mgb, Fe^Qis, 
Fe527(), Fes335 and ^"65406 as the indicators. Further, in order 
to improve the precision of the stellar properties using the 
minimization method we interpolated the tabulated indices of 
TMB03 on steps of 0.1 Gyr in age, 0.02 in metalUcity and 0.05 
in a/Fe. 

4.2. Ages, metallicities and a/Fe ratio profiles 

The model predicted lines strength profiles are shown in Fig. 
5 with solid lines. The derived quantities (age, metalhcity, ele- 
ment abundance, colors, mass-to-Ught ratios, and resulting x^) 
are shown in Fig. 7. Here we average together every four points 
to reduce the scatter. The theoretical lines strength indices match 
well the measured parameters in the inner regions of the galax- 
ies, although an acceptable value oix^ ~ 3 is obtained seldom 
(see Fig. 7), possibly indicating a global underestimation of our 
statistical errors, or, more probably, the presence of systematic 
effects. In particular, in the case of NGC1600, the large values 
of are mainly driven by the significantly low HyS (Fig. 5 and 
red solid fines in the right bottom panel of Fig. 7). The large val- 
ues ofx^ in the outer regions along the minor axis of NGC4125 
are caused by the lower Hp strength. In the case of NGC7619, 
the large values ofx^ in the outer parts along both the major and 



the minor axes are mainly due to large divergence of Fe50I5 
(green soUd fines). The SSP equivalent age, metaUicity and el- 
ement abundance ratio along the major and minor axis of our 
galaxies are also plotted on the left side of each panel. In each 
panel the solid lines indicate the fitting parameter profiles and 
the dotted lines show the 1 cr errors. As it can be seen in Fig. 7, 
all the galaxies are [a/Fe] overabundant, at the value of 0.25, 
0.1 and 0.3 dex for NGC1600, NGC4125 and NGC7619 respec- 
tively. The a/ Fe ratios are fairly constant in the central regions, 
with hints for a decrease in the outer parts. Table 7 fists the slope 
of Aindex/Ato^(r) of our galaxies which were derived by fitting 
a straight line to our measured data points. 

The nuclear (inside Re/8) parameters were derived by Trager 
et al. (2000b) and Denicolo et al. (2005). In order to compare 
to values available in the literature we calculate the mean pa- 
rameters inside R^/S. For NGC1600 and NGC7619 Trager et al. 
(2000b) derived the nuclear parameters [age, Z/H, a/Fe] with 
[8.6 ± 1.7 Gyr, 0.35 ± 0.05, 0.22 ± 0.02] and [14.8 ±2.3 Gyr, 0.2 
+ 0.03, 0.18 ± O.OIJ which compare to our value of [ 13.1 Gyr 
± 1.1, 0.35 ± 0.03, 0.317 ± 0.001] and [ 11.7 Gyr + 0.9, 0.25 
± 0.02, 0.284 ± 0.001] respectively. In the case of NGC4125, 
Denicolo et al. (2005) obtained [age, Z/H, a/Fe] with [5.9 ±3 
Gyr, 0.32 ± 0.1, 0.1] while we find [ 10.1 Gyr ± 0.6, 0.12 ± 
0.02, 0.095 ± 0.001]. 



4.3. Color and M/L profiles 

The Johnson broad band U-B, U-V, B-V, V-R, V-I, V-K, J-K, J- 
H, H-K color and M/L ratio in B, V, R, I, J, H, and K bands 
profiles were calculated using the Kjoupa initial mass func- 
tion Kroupa (1995) with the help of the SSP models (Maraston 
1998). For a clear presentation in the figures, we only show the 
U-B, B-V and V-I color and the mass to light ratios in B, V, I 
band in this papers. On the top right panels in Fig. 7, the blue 
and red solid lines stand for U-B and B-V color respectively; the 
black fines indicate the V-I color. For NGC4125 and NGC7619, 
the measured V-I color are also over plotted with solid dots. For 
NGC1600, unfortunately, we do not have reliable broad band 
color profiles and we take the aperture B-V color from Sandage 
(1973). As it can be seen from the plots, the models predicted 
color profiles agree reasonably weU with the measured colors ex- 
cept for small discrepancies. We note that Maraston et al. (2009) 
discuss the limitations of current SSP models in predicting accu- 
rate colors. In particular, the B-V colors of old (~ 12 Gyr), solar 
metalUcity SSP models are expected to be -0.09 mag bluer than 
what used here, while the differences are smaller for the redder 
colors. In addition, if our estimated ages are too small because 
of the presence of a second younger coinponent (Serra & Trager 
2007), then the predicted SSP colors would be biased and bluer 
than they should. Moreover, SSP models do not account for the 
presence of dust, which could be present especially at the centers 
of our galaxies. The middle right of Fig.7 shows the theoretical 
M/L ratio in B, V, I bands; the blue and green lines display the 
M/L ratios in B and V bands respectively; while the red lines 
present the I band M/L ratios. The red long dashed Unes show 
the dynamical M/L in the I band for NGC4125 and NGC7619, 
and in the R band for NGC1600 (see Sect. 5). The minimized 
of selected lines strength are presented in the bottom panel on 
right, red, blue and green lines present the minimize x^ of Hfi, 
Mgb and Fesois respectively; while the black lines show the total 
minimized All of the our giant elliptical galaxies show a red 
sharp peak, mainly due to metalUcity gradients (Peletier 1989). 
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Fig. 7. The best fitting SSP equivalent age, metallicity and element abundance ratio are showed in the left in each plot. The right top of each plot 
shows the Johnson broad band U-B, B-V, V-1 colors profiles, the blue and red solid lines stand for U-B and B-V color respectively; the black 
solid lines indicate the V-I color. For NGC4125 and NGC7619 the measured V-I color of Bender et al. (1989) are shown as dots. In the case of 
NGC1600, the aperture B-V color taken from (Sandage 1973) are indicated. The M/Lgv/ are shown in the mid of the right columns, the blue, green 
and red solid lines indicate the M/L in the B, V and I band respectively; the red dashed line shows the best-fit dynamical M/L in the I band for 
NGC4125 and NGC7619 and in the R band for NGC1600; the minimized of selected line strengths are presented in the bottom panel on the 
right, red, blue and green lines present the minimized of Hfi, Mgb and Fesois respectively; while the black lines show the total minimized 
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Fig. 7. Continued. 



Table 7. Line strength and total metallicity gradients: A index/ A log (r). 



Galaxy 


Range 


Mgb 


Fe5015 


Fe5270 


Fe5335 


< Fe> 


Fe5406 




Z/H 


a/Fe 


Age 


Name 




Slope 
+ 


Slope 
+ 


Slope 
+ 


Slope 

± 


Slope 

± 


Slope 

+ 


Slope 

+ 


Slope 
+ 


Slope 

± 


Slope 

± 


NGC1600 


0-0.5 


-0.787 


-0.339 


-0.179 


-0.196 


-0.187 


-0.207 


-0.125 


-0.342 


-0.042 


1.32 






0.067 


0.082 


0.036 


0.053 


0.038 


0.094 


0.028 


0.029 


0.012 


0.43 




0-1 


-0.742 


-0.537 


-0.268 


-0.222 


-0.245 


-0.236 


-0.172 


-0.288 


-0.044 


1.398 






0.069 


0.085 


0.037 


0.046 


0.036 


0.029 


0.033 


0.028 


0.009 


0.343 


NGC4125 


0-0.5 


-0.378 


-0.328 


-0.215 


-0.267 


-0.241 


-0.121 


0.224 


-0.225 


-0.037 


-3.527 






0.021 


0.036 


0.038 


0.020 


0.024 


0.017 


0.020 


0.013 


0.006 


1.234 




0-1 


-0.464 


-0.405 


-0.260 


-0.319 


-0.289 


-0.185 


0.173 


-0.235 


0.008 


-2.464 






0.023 


0.035 


0.033 


0.020 


0.021 


0.019 


0.020 


0.014 


0.006 


1.250 


NGC7619 


0-0.5 


-0.711 


-0.571 


-0.192 


-0.352 


-0.272 


-0.299 


0.078 


-0.343 


0.041 


-1.420 






0.049 


0.082 


0.035 


0.038 


0.028 


0.036 


0.021 


0.028 


0.010 


0.614 




0-1 


-0.787 


-0.953 


-0.280 


-0.437 


-0.359 


-0.264 


-0.026 


-0.327 


-0.034 


2.143 






0.048 


0.103 


0.035 


0.037 


0.028 


0.036 


0.039 


0.024 


0.011 


0.497 



5. Line strengths and the local escape velocity 

As discussed in the Introduction, we want to exploit the link 
that the galaxy formation process has established between the 
Mgb line strength and the local escape velocity to constrain 
the density profile of dark matter halos in the outer regions of 
galaxies. The Vesc is given by: 



(4) 



and therefore at each radius it is sensitive to the total mass den- 
sity profile up to very large radii; for example, in the case of a 
spherical density distribution p the gravitational potential <1) is: 



cD(r): 



M(< r) 



+ 47r pir'y'dr' 



(5) 



where M{< r) - 4:7r J^'^ r'^p(r')dr'. In the following we will cor- 
relate Mgb measured at the projected distance r with V^sc com- 
puted at the intrinsic distance r. As noted by Scott et al. (2009), 
we see no difference if instead we compute the projected quan- 
tity: 



v.. 



Jp.dz 



(6) 



where z is the line-of-sight. Moreover, having in mind the phys- 
ical scenario described in the Introduction, where more enrich- 
ment in Mgb is achieved if more gas is retained, one would ex- 
pect the difference AV between Vf „(r) and the average velocity 
of stars and gas at that location r to play the crucial role. In fact. 
Ay happens to be proportional to Vesc, because elliptical galax- 
ies are virialized systems. 
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5.1. Dynamical models 

We use the axisymmetric Schwarzschild's orbits superposition 
technique (Schwarzschild 1979) to derive the gravitational po- 
tential profiles. Thereby, also the stellar mass-to-Ught ratio, the 
internal orbital structure as well as the velocity anisotropy of the 
galaxies can be determined. Here, we only briefly present the 
basic steps, details about our implementation of the method are 
given in Thomas et al. (2004, 2005b). 

(1) : We first determine the luminosity density from the sur- 
face brightness profile. The photometric data for the modelling 
are deprojected into a three-dimensional, axisymmetric distribu- 
tion with specified inclination using the program of Magorrian 
(1999). 

(2) : The total mass distribution of the galaxies consists of the 
stellar mass density and a dark matter halo: 



p = Tv+pdm, 



(7) 



where Y is the stellar mass-to-Ught ratio, and v is the deprojected 
stellar luminosity density. 

(3) : The gravitational potential O can be derived by integrat- 
ing Poisson's equation once the total mass profile is obtained. 
Thousands of orbits are calculated in this fixed potential. 

(4) : The orbits are superposed to fit the observed LOSVDs, fol- 
lowing the luminosity density constraint. The maximum entropy 
technique of Richstone & Tremaine (1988) is used to fit the kine- 
matics data by maximizing the function: 



(8) 



where S is an approximation to the Boltzmann entropy and 
is the sum of the squared residuals to the kinematic data. The 
smoothing parameter a controls the influence of the entropy S 
on the orbital weights, see Thomas et al. (2004) for more details. 

Two different types of dark matter halo distribution were 
tried to recover the mass profiles. First, the NFW halo profile 
(Navarro et al. 1996): 



(r/r,)(l+r/r,)2' 



(9) 



where ps is the characteristic density of the halo and is a char- 
acteristic radius. One further defines the so-called concentration 
parameter c of the halo that is related to and the viral radius 
^"200 via Vs = r2oo/c. The potential generated by this density dis- 
tribution is given by: 



a)(r) = 4nGp,r,^^ Mr + -), 



(10) 



where G is the gravitational constant. The second halo used is 
the cored logarithmic (LOG) halo profile (Binney & Tremaine 
1987) , that reads as: 



pLOoir) 



(11) 



This gives a roughly constant circular velocity Vc and a flat den- 
sity core inside r < r^. The potential is given by: 



V2 

<D(r) = ^In(r2 + r2) 



(12) 



Both the NFW and the LOG potential give divergent mass pro- 
files when integrated to infinity. In our numerical implementa- 
tion the potential is computed up to 10 eff'ective radii and extrap- 
olated Keplerian at larger radii. In order to investigate further the 



physical extent of such halos we introduce a cut-ofl^ radius Vcut 
and modify the density profile given by Eq. 1 1 as follows: 

PcuAr) = pwoir) X exp - . (13) 



5.2. Model results for NGC1600, NGC4125 and NGC7619 

Photometric data: We produced isophotal fits separately on the 
HST and the ground based images using the code of Bender 
& Mollenhoff (1987). For NGC4125 and NGC7619, flie sur- 
face brightness profiles consist of HST F814W(I) filter images 
and SDSS i band images, scaled to I band. The photometric 
data of NGC4125 extended to 216 arcsec, about 3.5 times of 
effective radius. For NGC7619, the photometric data extended 
to 167 arcsec, more than 5 times of its effective radius. In the 
case of NGC1600, the photometric data came from HST images 
observed with the F555W(R) filter and scaled to the profile of 
Peletier et al. (1990). They extend to 200 arcsec, nearly 4.5 times 
of the effective radius. The photometric data are deprojected into 
the 3d luminosity distribution using the program of Magorrian 
(1999). Kinematic data set: The kinematics data both along the 
major axes and minor axes derived in Sect. 3.1 are used. 

We tested models at three different inclinations (i = 70, 80, 
90) for our galaxies, the best-fit models for all three galaxies 
are edge-on {i = 90). The dark halo parameters for the best- 
fit logarithmic halos are given in Tab. 8. They fit very well to 
the dark matter determinations of Coma early-type galaxies by 
Thomas et al. (2009). We also tested dark matter halos follow- 
ing the NFW halo profile, as in Thomas et al. (2005b). The 
best-fit NFW-halos for NGC1600, NGC4125 and NGC7619 re- 
quire concentration parameters c = 5, 5.5, 7.5 respectively. The 
mass-to-light ratios derived with NFW profile are similar to the 
mass-to-light ratios with logarithmic halo profile. Fig. 7 shows 
that the mass-to-Ught ratios derived from the SSP models with 
a Kroupa IMF are on average 45% lower than the dynamical 
ones, similarly to Cappellari et al. (2006). A better agreement 
is achieved if a Salpeter IMF is considered. The comparison 
with the mass profiles derived from X-ray measurements for 
NGC4125 and NGC7619 (Fukazawa et al. 2006) is good: our 
total masses within the last kinematic point are just 10% larger 
than the X-ray ones. For NGC1600 the discrepancy is larger, 
with X-ray masses underestimating the stellar kinematics ones 
by 30 (NFW potential) or 40 % (LOG potential). 

As a further step, we compute a number of models with vary- 
ing rcut from 5 to 10 times effective radii for each galaxy. The 
different cut-off radii of the dark matter halos do not have sig- 
nificant effects on the dynamical parameters (T, Vc, rc), nor the 
quaUty of the kinematic fit (see upper panel of Fig. 8, where the 
self-consistent model without dark matter is plotted at rcui - 
and the model with no cutoff at Vcm - 80 kpc), as soon as mod- 
els with rcut ^ 40 kpc are considered. However, they do have an 
effect, when we also take into account the correlation Mgh-Vesc- 

5.3. Mgb, the local escape velocity and the size of dark 
matter halos 

Next we explore the correlation between Mgb and the local es- 
cape velocity by fitting log Mgb and log Vesc to a straight-line, 
weighted by the errors. We computed errors on Vesc by consid- 
ering at every point the range of escape velocities allowed by all 
models giving ^ < ;rL,mi« + 1- According to eq. (5) tiie outer 
halo-profile has a considerable effect on Ve.u- and, hence, on the 
tightness of the correlation between log Mgb and log Vesc- Here, 
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Table 8. Halo circular velocity V^, halo core radius r^, average dark 
matter density < pdm > (within 2 r^fj ) and stellar mass-to-light ratio Y 
(last column: photometric band) for the best-fit dynamical models of 
NGC1600, NGC4125 and NGC7619 (no cutoff radius). 



Name 


Vc 


rc 


logio < Pdm > 


Y 


band 




(km/s) 


(kpc) 








NGC1600 


350+^'" 


Q1+52 


-2.1 ±0.5 


7.3 ±0.5 


R 


NGC4125 




32+7 


-2.2 ±0.1 


5.7 ±0.7 


1 


NGC7619 


600-00 


49+9 


-2.0 ± 0.2 


6.5 ± 0.5 


I 



we explore shallower halo profiles than probed by Scott et al. 
(2009) and use the correlation between Mg^ and Vesc to con- 
strain the size of dark matter halos. We define a new parameter 



Xtot , given by: 

2 _ 2 2 
Xtot ~ ^Msb ~^Xkin^ 



(14) 



where, xltgb quantifies the goodness of fitting log Mgb and log 

Vesc to a straight-line. Thexlf^ mentioned above, it quantifies 
the deviation between model and observed kinematics. 

We probed two different radial scalings: (1) correlating the 
measured Mg^(r) at radius r with the escape velocity Vesc(r) at 
the same radius r and (2) correlating Mgb with the escape ve- 
locity on the corresponding isophote (i.e. we compare the Mg^ 
measured at a distance b from the galaxy centre along the minor- 
axis with Vescib/ y/q) in the equatorial plane and the Mgb mea- 
sured at a distance a from the galaxy centre along the major-axis 
with Vescid V^' where q = 1 - e is the apparent flattening of 
the light-distribution). Scaling (2) gives a smaller and is dis- 
cussed in the following. Using the scaling (1), however, gives 
similar results with respect to the halo core radii. 

The lower panel in Fig. 8 shows xjot ^ function of the cut- 
off radius in kpc. The minimum scatter is achieved for rcut = 
60 kpc. This is 4.3 times the effective radii for NGC1600, 9.5 
times for NGC4125 and 8.5 times for NGC7619. Lower cutoff 
radii give slightly worse Mgb - Vg^c correlations, but the trend is 
significant only in combination with;^^^^.^. Cutoff radii larger than 
70 kpc have significantly higher As discussed above, this is 
driven only by the Mgb - Vesc correlation. 

The Mgb-Vesc correlation with the lowest is 

logMgZ; = (0.873 ± 0.008) x log V,,, - (2.105 ± 0.024) (15) 

(shown in the bottom part of Fig. 9 by the solid line), the corre- 
sponding ^^^^ is 592 (for 368 data points). The correlation for 
the best-fitting NFW halo reads 

logMgfo = (0.629 ± 0.004) x log Vesc - (1.314 ± 0.013) (16) 

(shown by the long dashed-dotted line), the corresponding 
is 1360. The regression obtained for models without dark matter 
is shifted to lower Vesc, but is steeper than reported by Scott et al. 
(2009). For comparison, in the upper plot of Fig. 9 we show 
the Mgb-Vesc relation obtained when the best-fitting logarithmic 
halo models without cutoff are used. Note that we get the same 
cut-off radii if we use [Z/H] - Vesc instead of Mgb-Vesc- 

Figure 9 indicates that the simple linear relation between 
log Vesc and log Mgb breaks down in the outer regions of the 
galaxies investigated here. There, the measured value of Mgb 
is much smaller than the linear relation determined in the in- 
ner parts of the galaxies would predict. Due to the relatively 
low number of points involved, the bending of the relation does 
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Fig. 8. The xli„ as a function of cut-off radius summed over all three 
galaxies is shown in the upper panel and the xl^gi, is presented in bottom 
panel. The dotted lines shows the;^'^ for the NFW halo profile, and the 
last point is the;^^ for the LOG dark halo profile without cutoff. 



not influence the result of Fig. 8: the same optimal cut-off ra- 
dius is found if we do not consider points with log Mgb < 0.5. 
However, this might be telling something about the formation 
process of ellipticals. It is reminiscent of several recent obser- 
vational and theoretical findings (van Dokkum et al. 2010, and 
references therein) that suggest that the outer steUar envelopes of 
bright ellipticals might be the result of later accretion of smaller 
objects on older central cores formed in violent collapsing pro- 
cesses at high redshift. In this case, the metallicity measured in 
the outer parts would be set be the (lower) escape velocity of the 
accreted objects, and not by the deeper potential well of the final 
giant elliptical. 

Figure 10 shows the total density profile derived for 
NGC1600 for the different cases discussed above (self- 
consistent, with a logarithmic dark matter halo with or without 
cut-off, or with a NFW halo). Clearly, the Vesc - Mgb correlation 
allows us to constrain the outer shape of the total matter den- 
sity, rather than delivering a precise determination of the radius 
where the dark matter halos terminate. In the next section we 
discuss this degeneracy in more detail. 



5.4. The degeneracy between outer halo slope and cut-off 
radius 

A steeper outer halo slope results in less mass in the outer re- 
gions and, thus, a lower escape velocity. Likewise, a lower cut- 
off radius reduces the outer mass and the escape velocity as well. 
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Fig. 9. The upper panel and bottom panel show the correlation between 
log(Mgb) and log(K(,,£.) derived using LOG dark matter halos without 
cut off and cut off at 60 kpc respectively. Different colour refer to dif- 
ferent galaxies. The full line is the corresponding best fit. The dashed- 
dotted lines is computed for NFW dark matter halos. The long-dashed 
line shows the regression obtained when models without dark matter 
halo are considered. The dotted line is taken from Scott et al. (2009) 
without dark matter halos. 



Consequently, a steeper density fall-off in the outer halo region 
could be partly compensated for by a larger cut-off radius. In 
the following we investigate this degeneracy quantitatively with 
spherical test models. 

The test models are set up as follows. We assume that 
the galaxy's light distribution can be described by a Hernquist 
sphere (Hernquist 1990). For simplicity we assume that the 
mass-to-light ratio equals one in appropriate units. Then, the 
stellar mass distribution p* reads 



Ma 
p. = — 



1 



In r r(r + ay 



(17) 



where M is the total stellar mass and a is the scaling radius 
(Hernquist 1990). In the following we use M = 3 x 10" M© and 
'"eif - 10 kpc (reft ~ 1.8fl for the Hernquist sphere). To investi- 
gate the degeneracy between slope and cut-off radius we seek a 
halo-profile whose slope can be adjusted conveniently. For ex- 
ample 



Pdm 



Po 



(r + ro)i 



(18) 



The four free parameters here are the normalisation po, a scal- 
ing radius ro (inside which the halo density eventually becomes 



Fig. 10. The total density profile derived for NGC1600 for the self- 
consistent model (dot-dashed line), the LOG dark matter halo cut at 60 
kpc (full line), the LOG dark matter halo without cutoff (dotted line), 
and the NFW dark matter halo (dashed line). 



constant), the cut-off radius rcut and the asymptotic logarithmic 
slope y (in the absence of a cut-off). For the rest of this section 
we fix two of these parameters according to empirical scaling 
relations for early-type galaxies. Firstly, the normahsation po is 
chosen such that 



P.(l-5refl^) =Pdm( 1-5 reft) 



(19) 



(Thomas et al. 2007b). Secondly, the halo scaUng radius is fixed 
to ro - 2.5 reft (Thomas et al. 2009). 

Since the goal is to investigate how much the outer slope 
is degenerate with the cut-off radius we first construct an input 
escape-velocity profile by fixing the remaining halo parameters 
■y and rcut to some values, e.g. y = 2 and rcut = 5 reft. We calcu- 
late the corresponding escape velocity profile Vesc from the total 
potential 



(20) 



where is given analytically for the Hernquist sphere 
(Hernquist 1990) and Odm is obtained by numerically integrat- 
ing equation (5) forpoM of equation (18). We project the escape 
velocity via Eq. 6. Note that the projection does not affect the 
results discussed below. For the aim to study the degeneracy be- 
tween halo slope and cut-off radius we assume that the projected 
escape velocity is given at 20 logarithmically spaced radii be- 
tween 0. 1 reft and 5 reft, with an uncertainty of A log Vescpr = 
0.05 (see Fig.9). 
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Fig. 11. Fitted cut-off radius r^ut against halo-slope. The input model is 
indicated by the red open circle. The filled black dots show the best-fit 
cut-ofi' radius for halo-slopes y = 2,2.4,2.8,3.2,3.6. Small grey dots 
indicate the region within Ax^ < 1 with respect to the input model. 
Left panels: input-slope y = 2 (isothermal-like); right panels: input- 
slope 7 = 3 (NFW-like). Top row: input r^ut = 5 r^ff; bottom row: input 
''cut = 12.5 r^jf . 



Can we fit the so constructed escape-velocity profile also if 
we make a wrong assumption upon the halo slope? To answer 
this question, the slope is reset, e.g. y = 3 and the remaining free 
parameter r^ut is fitted such that the;(f^ between the input escape- 
velocity profile and the corresponding profile with the new halo- 
slope is minimised. Fig. 1 1 shows that for a wide range of outer 
halo-slopes - we probed y e [2, 3.6] - the escape-velocity profile 
can be reproduced well (with ax^ < 1). As expected, a steeper 
halo-slope requires a larger cut-off radius to compensate for the 
lowered outer mass-density. Our first result is therefore that with 
the given uncertainties on the escape-velocity there is a large 
uncertainty in the outer halo-slope and cut-off radius, related to 
the degeneracy between the two. 

In our standard modeling, the halo-slope is y a: 2 (cf. Eq. 11). 
Next we ask how much this assumption can bias the 'measured' 
cut-off radii if the actual halo-slope is steeper. To this extent we 
repeated the above analysis for a wide range of input cut-off radii 
('"cut/'"efF 6 [5,20]) and an input slope of y = 3. The choice of 
y = 3 is driven by the fact that pure dark matter cosmological 
A^-body simulations predict a halo profile close to Eq. (9), i.e. 
y = 3. For each input model, we determined the 68% confidence 
region that would result under the assumption of y = 2. Fig. 12 
shows that cut-off radii reconstructed with y = 2 are always 
too small - as akeady discussed above. The figure also shows 
that if the true cut-off radius is increased by a factor of four, 
the reconstructed one (with y = 2) increases only by factor of 
two. In other words, if the steepness of the halo density profile 
is underestimated then (1) the cut-off radii are biased too low 
and (2) a relatively small uncertainty in the 'measured' rcut - 
obtained under the assumption of a y = 2 - may correspond to 



Fig. 12. Reconstructed halo cut-otf radii (assuming an outer halo-slope 
7 = 2) against true cut-off radii (the true halo-slope is 7 = 3). The points 
show the 68% confidence region of the reconstructed cut-off radii. The 
dotted line is the one-to-one relation. An underestimation of the steep- 
ness of the halo density profile results in an underestimation of the cut- 
off radius. 



a relatively large uncertainty in the true cut-oflF radii if the real 
value of y is significantly larger than two. 

6. Discussion and summary 

Accurate kinematic profiles extending out to 1.5 to 2 R,, along 
the major and minor axis for the three giant elliptical galax- 
ies NGC1600, NGC4125 and NGC7169 have been measured, 
together with 6 lines strength Lick indices. For NGC4125 we 
also detected gas emission along the major axis and measured 
its kinematics and equivalent width strength. From the compar- 
ison of the [NI]/l/l5197,5200/H/3vs [OIII] ^5007/1^8 diagnostic 
diagram to the previous work, we concluded that the emission 
region is probably caused by LINER-like emission. 

With the help of the SSP models, we derived the stellar pop- 
ulation parameters, the M/L ratios and the broad band colors of 
our galaxies. We found the galaxies NGC1600 and NGC7619 
have high metallicities, in contrast, the galaxy NGC4125 is just 
slightly above solar metallicity. The three objects have signifi- 
cant metallicity gradients. All of the galaxies are overabundant, 
NGC4125 just by 0.1 dex, and do not have significant gradients 
of element abundances along the axes. These galaxies have sharp 
red peaks at the center, reflecting the steep change of metallicity 
and in agreement with the models predicted color profiles. 

According to the simple element enrichment scenario, the 
a elements are mainly delivered by Type II supernovae explo- 
sions of massive progenitor stars, and a substantial fraction of 
Fe peak elements come from the delayed exploding Type la su- 
pernovae (Nomoto et al. 1984; Thielemann et al. 1996). Thus 
the a/Fe can be used as an indicator to constrain the forma- 
tion timescale of stars. Hence the absence of radial variations 
of a/Fe ratio in our galaxies likely suggests that there is no ra- 
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dial variation of star formation time scales. The radial metal- 
licity and lines strength gradients give one of the most strin- 
gent constraints on the galaxy formation. The galaxies that form 
monolithically have steeper gradients and the galaxies that un- 
dergo major mergers have shallower gradients. The mean metal- 
Ucity gradients for non-merger and merger galaxies derived by 
theoretical simulation in Kobayashi (2004) are a[Z/H]/A log(r) 
~ -0.30 ± 0.2 and -0.22 ± 0.2, respectively. The author found 
that the galaxies with gradients steeper than -0.35 are all non- 
major merger galaxies. As it can be seen from Table 7 the gra- 
dients A[Z/i/]/A log(r) of our galaxies are compatible with nu- 
merical simulations. At face value, there is a weak indication 
that NGC1600 and NGC7619 were formed through a monolithic 
collapse process, while NGC4125 was shaped via a (recent) ma- 
jor merger as also indicated by minor-axis rotation, the presence 
of gas, young central stellar ages, just mild a/Fe-overabundance 
and the unrelaxed appearance of the outer stellar envelope. 

Using the axisymmetric Schwarzschild's orbits superposi- 
tion technique (Schwarzschild 1979) with and without dark ha- 
los, we derived the local escape velocity of these galaxies. From 
the correlation between Mgb, and the local escape velocity, we 
further confirmed the suggestion of Franx & Illingworth (1990) 
that metallicity and lines strength are a function of local escape 
velocity. Moreover, we considered models with logarithmic dark 
matter halos of different sizes. The best fitting Mgb-Vesc relation 
comes from models which cut off the logarithmic halo at 60 kpc. 
Similar cut-off radii are obtained when using the correlation be- 
tween [Z/H] and Vesc but the scatter in this relation is larger, 
because it also depends on the stellar population age (see also 
Scott et al. 2009). We find that - at a given escape velocity - the 
youngest of our galaxies, NGC4125, is also the most metal-rich 
one, which could be explained within a merger scenario as star 
formation from enriched material. 

Our cut-ofF radii agree with the lensing analysis of cluster 
galaxies by Halkola et al. (2007) but not with the results from X- 
rays and strong-lensing (Humphrey & Buote 2010; Koopmans 
et al. 2009). Larger cut-off radii are compatible with the data if 
dark matter halos with steeper outer density profiles are consid- 
ered. 

We discover deviations from a linear correlation between 
log Mgb and log Vesc in the outer parts of our galaxies, where 
the measured Mgb are lower than predicted by the extrapolation 
from the correlation further inside. The physical link between 
Vesc and Mgb is established at the main time of star formation 
activity of the galaxy. Dry mergers happening after this episode 
will increase Vesc without increasing Mgb. Hydrodynamical cos- 
mological simulations by Naab et al. (2009) indicate that the 
outer parts of massive elUpticals (beyond ^ 5 kpc) might be 
dominated by stars which were bom in low-mass halos and ac- 
creted by frequent minor mergers. If this would be the dominant 
process for the formation of the outer parts of giant ellipticals 
one would expect such a break-down of the linear correlation 
between log Mgb and log Vesc as observed, because the outer 
Mg^ would come from stars which were formed in small halos 
with low Vesc (and accordingly low Mgb). 

We plan to further investigate the connection between line 
strengths and dark matter halos with a more extended galaxy 
sample. 
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